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We present a theory for the recombination of (charged) holons and doublons in one-dimensional organic Mott 
insulators, which is responsible for the decay of a photoexcited metallic state. Due to the charge-spin separation, 
the dominant mechanism for recombination at low density of charges involves a multi-phonon emission. We 
show that a reasonable coupling to phonons is sufficient to explain the fast recombination observed by pump- 
probe experiments in ET-F2TCNQ, whereby we can also account for the measured pressure dependence of the 
recombination rate. 
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Introduction - Femtosecond pump-probe spectroscopy is 
a powerful probe for the charge relaxation and thermalization 
phenomena in complex materials. These measurements can 
directly address and unveil the role of strong electron correla¬ 
tions, as well as the coupling to phonon degrees of freedom. 
Materials that behave as Mott insulators due to strong elec¬ 
tron Coulomb repulsion contain all the latter physics, and are 
therefore of high theoretical and experimental interest. It has 
been observed that photoinduced charges decay within the pi¬ 
cosecond range, i.e. well within the experimental resolution, 
but on the other hand orders of magnitude faster than in clean 
semiconductors with similar energy gaps. 

So far two classes of Mott insulators, investigated by pump- 
probe spectroscopy, revealed similar behavior. These are the 
layered undoped cuprates La2Cu04 and Nd2Cu04 HJ-O, and 
the quasi one-dimensional (ID) organic Mott insulators of the 
TCNQ family || 4 |, in particular ET-F2TCNQ iHJQ which will 
be the focus of our study. Both undoped cuprates and ET- 
F2TCNQ reveal ultrafast picosecond charge recombination 
with some similarities: (a) The charged carriers created by 
the pump pulse above the Mott-Hubbard (MH) gap are holons 
and doublons, and their recombination requires the distribu¬ 
tion of a large energy quantum (the MH gap A ^ leV) into 
several final excitations with smaller energy cq. At low den¬ 
sity of charges candidates for recipient bosons can be spin or 
phonon excitations, (b) The decay is exponential in time. This 
excludes bi- and higher-molecular processes involving inelas¬ 
tic collision of several ’free’ charge carriers, and implicitly 
reveals the existence of an intermediate bound state of a holon 
and a doublon (the MH exciton). In this respect a different ob¬ 
servation has been obtained on Ca2Cu03 from the ID cuprate 
family, which is known to have negligible excitonic effects m 
and thus shows a non-exponential decay || 9 |- 

From a theoretical viewpoint the challenge of understand¬ 
ing the charge recombination has analogies with the decay of 
the double-occupancy in ultracold bosons cni and fermions 
iniiii in optical lattices, where the decay rate E exhibits an 
exponential dependence on the ratio of the Coulomb repulsion 
U and the typical excitation’s energy scale eo- In the latter 
case the system can be described by a high-temperature state 
with a sufficient density of excited charges, so that the creation 
of particle-hole pairs in the compressible background m is 


the dominant decay channel, and eg is set by the kinetic en¬ 
ergy of recipient excitations, as observed also within DMFT 
03 . On the other hand, in real materials the final effective 
temperature is low, T U. For the case of 2 D undoped 

cuprates, which are antiferromagnets at low T, a theory has 
been presented M\l 5 \ where the fast charge recombination 
is explained via emission of spin excitations with the spin ex¬ 
change energy, eg ~ J, as the relevant excitation scale. Strong 
correlations and large J at the same time lead to a nontrivial 
origin of the s-type bound state of holon and doublon Il 6 l,i.e. 
the MH exciton, being the intermediate state essential for the 
exponential decay. 

In spite of similarities with 2 D Mott insulators, in quasi-ID 
Mott insulators the scenario involving spin-excitations cannot 
be effective neither for the MH exciton formation nor for the 
multi-boson emission due to the phenomenon of charge-spin 
separation. In the following we will show that a multi-phonon 
emission can be a viable recombination mechanism in ID or¬ 
ganic Mott insulators, somewhat specific to organic materials 
with energetic intra-molecular vibrations and strong electron- 
phonon coupling ifTTl [TSl . The mechanism bears similarity 
with recently proposed multi-phonon exciton decay in semi¬ 
conducting carbon nanotubes CSl while in standard semicon¬ 
ductors such a scenario seems to be inefficient ll^ . The pre¬ 
requisite is again the existence of the ID MH exciton II2TI . 
which can be stable in the case of longer range Coulomb re¬ 
pulsion. Since the photoexcited exciton is of the odd sym¬ 
metry we will show that its decay becomes allowed only due 
to the electron-phonon coupling. Finally we show that our 
scenario can explain the pressure dependence of the recombi¬ 
nation rate established recently for ET-F2TCNQ IjI. 

The problem is tackled as follows. By neglecting the re¬ 
combination term of the Hamiltonian we first compute the ex¬ 
citon, i.e., the lowest bound state in the sector with one dou¬ 
blon and one holon. We then use the Fermi’s golden rule in 
order to compute the decay of the exciton |T'i) into the man¬ 
ifold of states I'k™) that consists of the charge ground state 
with additional phonon excitations, which in the process of 
recombination receive the energy of the exciton. In principle 
virtual hoppings give rise to the spin exchange J = 
however, we shall neglected it in our calculation. Such ap¬ 
proximation is justified by the charge-spin separation specific 
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for ID, which makes the scattering of charges on spins inef¬ 
fective (since the hopping of holons and doublons only shifts 
the spin background), and by the hierarchy of energy scales 
J ujo,t,V,U for typical organic materials. 

If the Mott gap A is of the order of several phonon fre¬ 
quencies ujQ, A « nuJo, the n-phonon contribution deter¬ 
mines the matrix element in the Fermi’s golden rule expres¬ 
sion. The electron phonon-problem is controlled by two di¬ 
mensionless parameters, the coupling strength ^ = A^/wg 
and the adiabaticity t/ujQ. In the most general case, com¬ 
puting the exciton in the presence of electron-phonon inter¬ 
action is not possible analytically. To generate admixture of 
n ^ 1 phonons to the exciton, the coupling strength ^ must 
be treated to higher orders. On the other hand, at least the 
limit t/ujQ <C 1 is a valid starting point for molecular vibra¬ 
tions in organic crystals. In this case it is convenient to rewrite 
the Hamiltonian using a unitary Lang-Firsov transformation 
e^, which measures the phonon coordinate with respect to the 
equilibrium position for a given charge configuration. If the 
transformed exciton state |^i) is expanded in phonon number 
states, |T'i) = ..., the zero-phonon 

state is already the leading contribution in t/uJo so that 

additional phonon-dressing can be neglected. Using this ap¬ 
proximation we will derive a compact expression for the re¬ 
combination rate F, 



which can easily be compared with the experiments, taking 
the hopping t and the nearest neighbor interaction V from in¬ 
dependent measurements. 

The Model - As a model for the charge recombination in 
organic Mott insulators we consider the ID extended Hubbard 
model, where in addition to the local Hubbard repulsion U 
and the nearest-neighbor electron hopping, a nearest-neighbor 
Coulomb repulsion U > 0 is included. The latter is essential 
to stabilize the exciton state in ID ED- The Hamiltonian is 
split in the hopping Ht of doublons and holons, the recombi¬ 
nation term Hrc and the interaction term Hjj, which are writ¬ 
ten as 


Ht = 

Y (4,sd3,s 

- dlsh],s+Hc.), 

(2) 


(i3),s 



Hrc — 


s “1“ “1“ H-C.), 

(3) 


{i3),s 



Hu = 

TT^ nf-f nf 

^ 2 

+ V Y'dihj, 

(4) 


(ij) 


with holon and doublon creation operators hi ^ = Ci,s(l — 
= clgTii^s, holon and doublon density operators 

;, and hi = nf — nf. 


-iV /?'/?■ — AV 

“ 2 "-i — 2 2.^1 


dl „di 


Here {ij) denotes nearest neighbor pairs and s the spin oppo¬ 
site to s. In addition a generally nonlocal coupling between 


the charge density and dispersive phonons is introduced. 

Hep — ^ ^ AgC (flg ~\- (l—q)lflj^ Hpli — ^ ^ UJqO^^dq . (5) 

3,1 9 

Lang-Firsov transformation- The derivation of the stan¬ 
dard Lang-Firsov transformation for the present case follows 
Ref ifTSll and is presented in the Supplemental Material. The 
exact transformed Hamiltonian is given by 

iTo = -i (dlsdjse^^^'e^^^ - -f H.c.) 

__ _ _ __ 

3 (ij) 9 

Tfrc ~ f ^ ^ ^hisdjs^ ^*6 -f hjgdisG C -f H.C.^ , 
(i3),s 

(7) 


where A^-, = ~ ^ is a phonon 

creation term, and U, V, and t are renormalized interac¬ 
tion and hopping parameters. Corrections to the bare U, V 
are given by [/ = U — cq and V = V — ei with Ci-j = 
2 cos{q{i — j)), while longer range interaction 

shall be neglected. Below we will express all results in terms 
of the renormalized parameters U, V, and i, which are deter¬ 
mined experimentally by a fit to the linear absorption spec¬ 
trum. 

Exciton ground state - We now construct the ground state 
for the Hamiltonian Hq. To neglect additional phonon dress¬ 
ing as explained above, Hq is projected to the phonon vacuum, 
= \Oph){Oph\Ho\Oph){Oph\- We first construct a basis of 
all holon-doublon states with an arbitrary spin configuration 
of the remaining sites, analogous to the squeezed spin state 
Ea. For a given spin configuration a = {ai, ...,aL- 2 } (with 
aj =t, i), we define ) as the state obtained by distribut¬ 
ing the spins a on lattice sites {1,..., L} \ For m < j 


=c^ 


l,fTi 


■ • c. 


m—l,cr^_i 


J 

S + l.o-j-l 




A 

JO), 




( 8 ) 


and for j < m analogous. We then define the holon-doublon 
state jo;™) by placing a doublon at site j, 


k.r) 


0 if m = j. 


(9) 


These are used to define the state with a holon-doublon pair 
1$™) and the state with two holons |<i)™’-' ) on an arbitrary spin 
background which is a superposition or mixture of configura¬ 
tions a, 


\^f) = Y = E (10) 

(T G 


One can see that the Hamiltonian H^^ does not mix different 
configurations ct, because nearest neighbor hopping of a holon 
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or doublon implies a shift of the spin background, which is 
implicit in the definition (|^ for Ip:}) |p:}=^^) or (jp;}) 

The action of the Hamiltonian is thus obtained by 

(Hf -1/ + = -i ^ (!<!>;+“) - |<i>;.+j) 

a— ±1 

( 11 ) 

To determine the ground state we start from a partial Fourier 
transform with respect to the average position, 

^ j 

With this the action of the Hamiltonian becomes 

(#(0) -u + = -2ii, Y. (13) 

Q=± 1 


Similar, PoHrc^YKh) = -i e4.+i|0p^). 

In summary, 

X (e^i+i.j - (18) 

When this is inserted into Eq. recombination rate is ex¬ 
pressed as 

r = P/32^Pdrf(A), (19) 

d 

with a spin structure factor 

9d = (20) 

j 


where tq = fsin(p/2). There is a continuum of states in the 
energy window E G \U — Aiq,U + Atq]. For parity- 

even and odd bound states are degenerate. We can restrict the 
analysis to the odd states, which can be created by the optical 
dipolar transition, and thus make the ansatz 

1^®) = Pi = (14) 

l>0 

The ground state is found for q = n with Pi = Po{2t/Vy 
and El = U — V — 4f^/y, which lies below the continuum 
for V > 2t. Without the electron magnon coupling, the exci- 
ton is decoupled from the spin background, i.e., excitons for 
different spin wave functions are degenerate. 

Exciton decay - Similarly to the problem of exciton decay 
in 2D im [T5| we establish the recombination rate using the 
Fermi’s golden rule 

r = 27r^ - El), (15) 

ra 

for transitions from previously determined exciton 
Eq. (tg, into the charge ground state with additional phonon 
excitations |4'™) via the recombination operator Hrc, Eq. (|^. 
If written in an integral form ifTSll . Eq. ([T5|) becomes 


and a phonon emission factor 

nOO 

= 2Re / (0p^|(e^.'+i,.' -e^.'-'+i) x 

X - e^Ij+i)|0p;i). (21) 

Spin structure factor - From Eqs. ( |20l i and one can 

see that go — 2 and gi = p-i = 1 for an arbitrary spin 
configuration. For d > 2, Eq. ® implies that for any spin 
configuration a, {—l)'^{apd+i-^^g.^^^gp+d,j+d+i'^ equals 1 
if the spins (aj, ...,aj+d-i) form an antiferromagnetic se¬ 
quence ■•■) or ■■■)’ 0 else. We thus 

have gd^o = 1 for a perfect Neel antiferromagnet, and gd = 0 
for |d| > 2 for a spin-polarized background. For a general 
finite temperature state we expect an exponential decay of the 
correlations with distance. 

Boson emission factor - The matrix element in the boson 
factor ( |2 T] i can be evaluated straightforwardly, which is done 
in the Supplemental Material. We obtain 

poo 

Tf{A)=8Re drP^^ X (22) 

Jo 

X sinh(2y^ cos(dq)(l — cos. 

q Q 


r = 


poo 

2Re{¥y'>\HrcPo (16) 

Jo 


where A is the charge gap, Pq is the projection to the zero 
charge sector, and Hph is the only part of Hq which is active 
in the zero charge sector. 

We first evaluate the application of Hrc on the exciton. 
Starting from the expression Eq. (14l, we can restrict 

the application of Hrc to nearest neighbor terms 1 = 1, 


PoHrc\^+i)\0ph) = -t 5,1$^’^+') e^.+i-lOp/,), (17) 






where Sj = creates a spin singlet on 

sites j,j + 1, previously occupied by a holon-doublon pair. 


The argument of the sinh may be written in a convenient way 
as an integral 2 J duje~^‘^'^ fd(oj), with the boson coupling 
function 

fd{uj) = E cos(d 9 )(l - cosq)6{uj - ujq). (23) 

9 ^^9 


The zeroth and first moments gd = f dujfd{uj), UdVd = 
J dujLofdiijj) of these functions are related to the phonon- 
mediated long-range interaction parameters id via 2U^d = 
[id — \^d-i — ^£d+i)- The time-integration in Eq. (22 1 can 
be performed numerically for any kind of dispersions A 


q,^q, 


however, for a fixed function fd and A/ujq oo, one can 
use an argument related to the central limit theorem to show 






4 




Figure 1: (Color online) Comparison of boson emission factors 
for different d, obtained by numerical integration of Eq. \22\ for dis¬ 
persions ujq = ijJo + 5u! cos{q), Xq = Xl^fL (dashed lines) and from 
Eq. l |24| l (solid lines). Parameters wo = 1, do; = O.Icuq, ^ = 0.3 are 
used. 


Figure 2: (Color online) Recombination time as a function of pres¬ 
sure p calculated from the Eq. l[^ using the experimental parameters 
wo = 0.23eV,U = 0.845eV, and f(p) G [0, 04,0.OBjeF, 
V{p) G [0.12,0.16]eV for coupling strengths ^ = 0.25,0.27, 0.3. 


that in the lowest order the result depends only on the zeroth 
moment rjd, as presented in the Supplemental Material. The 


integral (22 1 can then be approximated with 


rf (A) = 


M / 27 r / A y 

Tjd V Awo V2e|pd|wo^ 


(24) 


where wg is the typical phonon frequency. Expression (24i is 


obtained also by the saddle point approximation for a Gaus¬ 
sian /d(w) with the same zeroth moment rjd ifTSl . 

Typically, phonons are weakly dispersive and electron 
phonon interaction not long-ranged, so that |eo| |ei| 
|e 2 |--- Then it suffices to take into account only the |d| < 1 
contributions, as demonstrated in Fig. [T] for dispersions Xq = 
a/Vt, Wg = Wg + 5uj cos{q), showing results of numerical 
integration of Eq. ( [22] l (dashed lines). To boost the conver- 

' ’ _ • _ • _ 2 2 /r) 

gence an additional smoothening e e 

with T] = 0.2a;o has been used in Eq. ( |22] i, which can phys¬ 
ically correspond to higher dimensionality of phonons. The 
final expression for the recombination rate, Eq. ([T]), which is 
relevant for the comparison with experiments, is thus obtained 
by restricting Eq. ([T^ to the |d| = 0,1 contributions with spin 
structure factor go — 2,g±i = 1, and using approximations 
fig « fli « wg, 2wg?7g « Cg, 2uJogi « —eg/2inEq. (241 with 
eg « 2a;g^ expressed via coupling strength ^ De¬ 

pendence r2^(A), Eq. (24 1 , with pd approximated as above is 
for relevant terms d = 07l shown in Fig.(solid lines), dis¬ 
playing agreement with numerical integration. The prefactor 
(1/2 — 2P(V"^) in Eq. ([^ comes from Pf. 

Comparison with experiment - Finally we compye the 
recombination times obtained from Eq. (ID with 

experimentally measured ones Q. All quantities but the 
strength of charge-phonon coupling are set by the exper¬ 
imental data: wg = 0.23ey 1231 . U = 0.845ey, while 
i{p) G [0,04,0.06]eE, V{p) G [0.12,0.16]eE are specified 
functions of pressure ||71, and A « i/-E-4P/'^-Fig.§dis- 
plays Tr as a function of pressure p for three different values 
^ = 0.25,0.27,0.30, showing that p « 0.27 is consistent with 
the experimentally measured recombination times ||7], yield¬ 
ing the electron-phonon coupling A = UJo^/p = 0.12eE. The 


latter has been measured and calculated for a similar organic 
material, finding A G [0.05eE, O.leE] ifTTll . confirming that 
the electron-phonon coupling needed to reproduce the experi¬ 
mental results is indeed realistic. 

Conclusions and discussion - The central result of our 
study is that fast charge recombination observed recently in 
quasi-ID organic Mott insulators ||7l can be explained via cre¬ 
ation of phonon excitations, which can be for the material con¬ 
sidered (ET-F 2 TCNQ) identified as molecular vibrations. Due 
to the charge-spin separation in ID systems and hierarchy of 
energies J < wg in materials addressed, spin excitations are in 
contrast to 2D systems an inefficient decay channel and were 
neglected in our analysis by setting J —0. Motivated by 
the experimentally observed exponential decay of charge den¬ 
sity we derive the recombination rate based on the assumption 
that holon and doublon initially form a bound state - exciton, 
which is odd under the parity transformation (therefore opti¬ 
cally accessible). Still, the transition into the charge g.s. with 
even parity is allowed due to the coupling to phonons. We 
established the charge recombination rate using the Fermi’s 
golden rule, showing approximately exponential suppression 
with the number of phonons emitted in the process; an ob¬ 
servation common to several doublon decay processes M- 
MM with different recipients of doublon energy. 

The experimentally established frequency of the relevant 
vibrations ll23l is much larger than the one of typical lat¬ 
tice phonons, making the recombination mechanism some¬ 
what specific for organic insulators. To explain recent ex¬ 
periments on ID cuprates (Ca 2 Cu 03 ) ii with smaller typi¬ 
cal phonon frequencies and negligible excitonic effect some 
modification of the mechanism might be needed and remains 
as a future challenge. One should note that to assist a 
proper dissipation of energy in the case considered the vibra¬ 
tions must be at least partially dispersive or coupled to other 
modes. Even though our derivation focuses on ID phonons, 
it is straightforward to generalize it to more realistic three- 
dimensional electron-phonon coupling. Recognizing the role 
of the electron-phonon coupling in the recombination mech¬ 
anism we see the recombination measurements as an indirect 
way to establish its typically elusive strength at least for this 













5 


Okamoto. This work has been supported by the Program 
PI-0044 and the project J1-4244 of the Slovenian Research 
Agency (ARRS). Z. L. is supported also be the L’Oreal - UN¬ 
ESCO national scholarship ’’For Women in Science”. 

The authors acknowledge helpful discussions about the 
experimental result with S. Kaiser, M. Mitrano, and H. 


class of materials. 


Acknowledgments 


[1] K. Matsuda, I. Hirabayashi, K. Kawamoto, T. Nabatame, T. Tok- 
izaki, and A. Nakamura, Phys. Rev. B 50, 4097 (1994). 

[2] H. Okamoto, T. Miyagoe, K. Kobayashi, H. Uemura, H. Nish- 
ioka, H. Matsuzaki, A. Sawa, and Y. Tokura, Phys. Rev. B 82, 
060513 (2010). 

[3] H. Okamoto, T. Miyagoe, K. Kobayashi, H. Uemura, H. Nish- 
ioka, H. Matsuzaki, A. Sawa, and Y. Tokura, Phys. Rev. B 83, 
125102 (2011). 

[4] H. Uemura, H. Matsuzaki, Y. Takahashi, T. Hasegawa and 
H. Okamoto, Journal of the Physical Society of Japan 77, 113714 
(2008). 

[5] H. Okamoto, H. Matsuzaki, T. Wakabayashi, Y. Takahashi and 
T. Hasegawa, Phys. Rev. Lett. 98, 037401 (2007). 

[6] S. Wall elai, Nat. Phys. 7, 114 (2011). 

[7] M. Mitrano et al, Phys. Rev. Lett. 112, 117801 (2014). 

[8] M. Ono, K. Miura, A. Maeda, H. Matsuzaki, H. Kishida, Y. 
Taguchi, Y. Tokura, M. Yamashita, and H. Okamoto, Phys. Rev. 
B70, 085101 (2004). 

[9] H. Matsuzaki, H. Nishioka, H. Uemura, A. Sawa, S. Sota, T. 
Tohyama, and H. Okamoto, Phys. Rev. B 91, 081114 (2015). 

[10] A. L. Chudnovskiy, D. M. Gangardt and A. Kamenev, Phys. 
Rev. Lett. 108, 085302 (2012). 

[11] N. Strohmaier, D. Greif, R. Jbrdens, L. Tarruell, H. Moritz, 
T. Esslinger, R. Sensarma, D. Pekker, E. Altman, and E .Dernier, 


Phys. Rev. Lett. 104, 080401 (2010). 

[12] R. Sensarma, D. Pekker, E. Altman, E. Dernier, N. Strohmaier, 
D. Greif, R. Jordens, L. Tarruell, H. Moritz, and T. Esslinger, 
Phys. Rev. B 82, 224302 (2010). 

[13] M. Eckstein and P. Werner, Phys. Rev. B 84, 035122 (2011). 

[14] Z. Lenarcic and P. Prelovsek, Phys. Rev. Lett. Ill, 016401 
(2013). 

[15] Z. Lenarcic and P. Prelovsek, Phys. Rev. B 90, 235136 (2014). 

[16] T. Tohyama, Phys. Rev. B 70, 174517 (2004). 

[17] A. Girlando, The Journal of Physical Chemistry C 115, 19371 

( 2011 ). 

[18] H. Matsueda, S. Sota, T. Tohyama and S. Maekawa, Journal of 
the Physical Society of Japan 81, 013701 (2012). 

[19] V. Perebeinos and P. Avouris, Phys. Rev. Lett. 101, 057401 
(2008). 

[20] P. Yu and M. Cardon, Fundamentals of semiconductors: 
physics and materials properties (Springer, Berlin, 19996). 

[21] F. H. L. Essler, E. Gebhard and E. Jeckelmann, Phys. Rev. B 
64, 125119 (2001). 

[22] M. Ogata and H. Shiba, Phys. Rev. B 41, 2326 (1990). 

[23] S. Kaiser effl/., Sci. Rep. 4, 3823 (2014). 

[24] R. Sensarma, D. Pekker, A. M. Rey, M. D. Lukin and E. Dern¬ 
ier, Phys. Rev. Lett. 107, 145303 (2011). 



Supplemental Material 
to 

Exciton Recombination in One-Dimensional Organic Mott Insulators 


6 


Zala Lenarcic^, Martin Eckstein^ and Peter Prelovsek^ 

^J. Stefan Institute, SI-1000 Ljubljana, Slovenia 

^Max Planck Research Department for Structural Dynamics, University of Hamburg-CFEL, 22761 Hamburg, Germany 
^ Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia 


I. LANG-FIRSOV TRANSFORMATION 

The unitary Lang-Firsov transformation is obtained by 
choosing the operator S of the form 

S = '^al{al-a-q)nj, (25) 


are phonon-induced long-range interaction parameters. For 
the transformation of the hopping we introduce a different rep¬ 
resentation of Eq. ( |25] l, 

S = ^Pjnj, (34) 

3 


3,q 


i.e., a coherent displacement of the phonon coordinate de¬ 
pending on the holon-doublon configuration. The parame¬ 
ters are chosen such that the direct coupling term TJep in 
the transformed Hamiltonian H = e~^He^ is eliminated, 
i.e., Hep — [S', Hpfj] = 0, which is achieved by the choice 
jojq. To see this we first write 




XoNo 


U3„ 


% + 


LOn 


(26) 


with Nq = = N^_q. In this representation, the 

terms for individual momenta commute, so that one can easily 
compute the transformation of Hph and Hep. Using expres¬ 
sions for a coherent state shift of bosonic operators. 


P3 =-p\=H{-^ ^“ 9 ) ■ (35) 

q \ iOq OJq V 

It is important to note that the operators ftjPj andp/n; do com¬ 
mute, because the density operators commute, and we have 

[pj,pi] sin {q{l - j)), (36) 

9 ^^9 

which vanishes due to Xq = Al^, ujq = LO-q. With this we 
have (because fijdjs = 0, djsfij = djs) 

e-^djse^ = djs . (37) 


we have 


gaat-a-a^g-aat+a*a 

gaat-a*a ^tg-aat+a*a = _ q,* ^ 


e 


e 







XqNq 

LOq 

X-qN q 

LOq 


An analogous calculation for the holon operator gives. 


(27) 




(28) 


e ^hjsG^ =hjs€ . 

(38) 


Equations ( 

7T\ and (|38]l can now be used to transform the hop- 


ping part (|2 

Jof the Hamiltonian, 


(29) 






e~^h^jehj>se^ = hJj^hj>sePH~PU, 

(39) 

(30) 


e~^d}jedjise^ = dj,,dj'se~P'’ePP. 

(40) 


Hence, writing Eq. Q as Hep = J2q {XqalNq-i-X_qaqN_q), 
we have 

e~^{Hph + Hep)e^ = 

= y^Wqa\aq-y^^^^NqN_q ( 31 ) 

q q ^ 

= ^ > (32) 

9 33' 

where 

^3-3' =255 cos {q{j - /)) = ( 33 ) 

- LOq 

n ^ 


In these expressions, electron-phonon interaction is present 
through the factors and e~PHPt'. 

Below we will study the action of the Hamiltonian Hq, in 
particular within the zero-phonon sector. For this purpose, 
it is convenient to rewrite all phonon-operators, especially 
the terms ePH~Pt' and e~PHP3' in Eqs. (|4gi and in a 
normal-ordered form, so that they give zero when acting on 
the phonon vacuum. With Eq. we define 

p^^A]-A„ A] = -Y,e-^‘^^^al. (41) 

^^ UJa ^ 

q ^ 

Using the Baker Hausdorff relation 
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and one then gets 

ePjg-Pj' = g^jj'^ (42) 

g-PjgPj' _ g-^jj'g^jj'g-Cj-j'^ (43^ 

(44) 

e,-,- = - l[Al,A,] - l[A]„A,,] (45) 

= Y1 “ cos(g(j - j'))]. (46) 

q ^<1 


The exponentials can be evaluated using Eqs. ED and ( |44l i, 
yielding 


4 sinh (2 


|A, 


9 




cos(dg)(l — 0039)6 *“ 9 '^). (53) 


Now we argue how the time integration in Eq. (22 1 can be 
approximately evaluated. While the precise form of the bo¬ 
son coupling function /^(a;) is often not known, for optical 
phonons it is centered around some frequency ijJq, and could 
be approximated with a Gaussian form 


When the hopping Hamiltonian is projected to the phonon 
vacuum, only a renormalization of the hopping by a factor 
(the Eranck-Condon factor) remains. Using Eqs. ( [39l l, 
( |40l i, ED’ ED’ ED transformed hopping Hamilto¬ 
nian Ht = e~^ Hte^in the Lang-Firsov representation is 
obtained, and together with the interaction term Hint reads 
Hq = Ht + Hint, 




j g ^ii 




Hint = E Y E ^ E 

9 


2 ^ 


H.c. 

(47) 

(48) 


where U = U — eg and V) = i5/,iU — e/ and 

i = (49) 


are renormalized interaction parameters. The parameters i 
and Vi are determined experimentally, by characterization of 
the exciton in linear absorption. For simplicity, we thus in¬ 
clude only the nearest neighbor interaction Vi = V in the 
simulation, other parameters are small and experimentally not 
known, while the calculation can be straightforwardly ex¬ 
tended to longer range hopping. 

Similar to Eqs. ( |39| l and ( |40l l, the recombination term Hrc 
reads 


Hrff, t 




hisdjgC ’96 hjgdisC 9’6 -f H.c. 


(50) 


II. EVALUATION OF THE BOSON FACTOR 

First we evaluate the integrand of Eq. ( |2T] ). It can be written 
as (for d = j - j') 

{Qph\{e^^'AA - - e^^.’‘+0|0p,.), (51) 

where the time argument of Ajj: is due to the evolution with 
Hph, which simply amount to replacing aq by ag6“*“9'^. Us¬ 
ing the Baker-Hausdorff relation 6^6^ = 6^6^61^’^! we can 
normal-order the bosonic operators, which gives 


5d(w) = 


Vd 




(54) 


V^CTd 

for which the integra l ([22] l can be established using saddle 
point approximation ISIL yielding decay rate of form ( |24) l. 
For a more general coupling function centered at a typical fre¬ 
quency Wo an argument related to the central limit theorem can 
be used to show that in the limit A/wq —> 00 the decay rate 
is up to the leading order determined by its zeroth moment 
Vd = J dujfdiuj) as written in Eq. (24l. 

We define 


/ OO 

duje~ 

-OO 




7d(^), 

f —OO 

/d(w) = —fd{oJ + OJo), 

Vd 

|CJ 5 


(55) 

(56) 

(57) 


where * • • • * is the m-fold convolution of fd{^)- The expres¬ 
sion for the recombination rate Eq. (|22]l can then be simplified 


T^’^{A)=8Re J dre*^’’sinh ^2 J duje~^‘^'^ fd{(jj)j 


= 4 E 


dTe^^^ifd{T)y 


(58) 

(59) 


I odd 


= 8^ M y y-)(A - mwo). (60) 

Vd ml 

m odd 

The TO-fold contribution R^\x) of the normalized func¬ 
tions becomes broadly centered around x = 0 for large m. 
The dominant contributions to the sum will come from terms 
around n = A/wq, therefore we can make further approxima¬ 
tions 


rf(A) 


Stt M (2^ y y»)(A - mwo). (61) 
m odd 


where n = [nj stands for the integer part of n. Using that for 

(n) 

large n with broad 

E ~ ^ / dujT^J'\A - w) (62) 


m odd 


1 

2a;o 


( 63 ) 
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and the Stirling approximation n\ « y/2nn{nf e)'^, we finally 
obtain the compact expression ([24)i, 


rf(A) 


\vd\ (2|pd|) 


Stt _ 

Vd n\ 

2wo 

Stt \r]d\ 1 Wo 

f 

2wo r]d V 27rA 

VectJQ/ 


A 


(64) 

(2|77d|)^ (65) 

^ ( 66 ) 
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